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Abstract 

We report results of the 2006 April multi- wavelengths campaign of SS 433, focusing on X-ray data 
observed with Suzaku at two orbital phases (in- and out-of- eclipse) and simultaneous optical spectroscopic 
observations. By analyzing the Fcxxv Ka lines originating from the jets, we detect rapid variability of 
the Doppler shifts, dz/dt w 0.019/0.33 day~^, which is larger than those expected from the precession 
and/or nodding motion. This phenomenon probably corresponding to "jitter" motions observed for the 
first time in X-rays, for which significant variability both in the jet angle and intrinsic speed is required. 
From the time lag of optical Doppler curves from those of X-rays, we estimate the distance of the optical 
jets from the base to be w (3 — 4) x 10^^ cm. Based on the radiatively cooling jet model, we determine 
the innermost temperature of the jets to be Tq = 13±2 keV and 16±3 keV (the average of the blue and 
red jets) for the out-of-eclipse and in-eclipse phase, respectively, from the line intensity ratio of Fcxxv Ka 
and Fcxxvi Ka. While the broad band continuum spectra over the 5-40 keV band in eclipse is consistent 
with a multi-temperature bremsstrahlung emission expected from the jets, and its reflection component 
from cold matter, the out-of-eclipse spectrum is harder than the jet emission with the base temperature 
determined above, implying the presence of an additional hard component. 
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1. Introduction 

SS 433 is a unique Galactic binary system that exhibits 
continuous relativistic jets, and hence is an ideal tar- 
get to study fundamental problems of astrophysical jets, 
including the acceleration/collimation mechanism, struc- 
ture, and influence on the environment (for a review, re- 
fer to e.g., Margon 1984; Fabrika 2004). The intrinsic 
speed of the jets is 0.26c (c is the light speed), and the 
jet axis precesses with a period of 162 days. The optical 
and X-ray spectra show several pairs of Doppler shifted 
lines from the bipolar jets. The strong emission lines, such 
as Ha and lines from ionized heavy elements, give direct 
evidence that the jets mainly contain baryonic plasmas. 
Most recently, Kubota et al. (2010) constrains the mass 
of the compact object to be 1.9-4.9 Mq, indicating that 
the supercritical accretion star is most likely a low mass 
black hole, although a possibility of a massive neutron star 
cannot be ruled out. The inclination angle and orbital 
period are accurately measured to be i = 78.8° (Margon 



& Anderson 1989) and P=13.082 days (Goranskii et al. 
1998), respectively. 

Soon after the discovery of SS 433, many efforts were 
done to determine the kinematic model and parameters 
of the jets, by tracing the Doppler shifts of emission lines 
in the optical spectra. Margon et al. (1979) for the first 
time proposed a kinetic model, where the jets of about 
1/4 the light speed make precession, with a period of 
162.15 days, according to the latest value obtained by 
Gies et al. (2002). In superposition to the "precession" 
motion, another regular periodicity of the Doppler shifts 
is observed as a "nodding" motions with an amplitude of 
Az ~ 0.01 and a period of 6.28 days. They are caused by 
the tidal torques from the companion star. Furthermore, 
faster random variability of Doppler shifts were some- 
times reported (see e.g., Ciatti et al. 1981; Katz & Piran 
1982; lijima 1993; Collins & Garasi 1994). These phe- 
nomena are called as "jitter" , which contains an impor- 
tant clue to understand the collimation mechanism of the 
jets. The origins of "jitter" have not been understood. 
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in particular whether they are caused by the change of 
the speed or orientation due to insufficient observational 
quality (CoUins & Garasi 1994). 

SS 433 is also a target of great interest in high energy 
astrophysics, and has been observed by many satellites 
from early on (e.g., with EXOSAT, Watson et al. 1986; 
Tenma, Matsuoka et al. 1986; Ginga, Kawai et al. 1989). 
The ASCA satellite, which carried X-ray CCD cameras 
for the first time, detected many pairs of Doppler shifted 
emission lines from ionized metals such as Si, S, Ar, Ca, Fe 
and Ni originating from the twin jets (Kotani et al. 1994). 
By detailed analysis of the line intensities of Fexxvi Ka 
to Fexxv Ka, Kotani et al. (1996) estimated the length 
of the X-ray jets and temperature of the jet base to be 
~ 10^^ cm and ~20 keV, respectively. The emission lines 
are resolved with the Chandra HETGS, which are found 
to have Doppler widths of wlOOO-5000 km s"^ (Marshah 
et al. 2002; Namiki et al. 2003; Lopez et al. 2006). The 
broad band continuum (up to 100 keV) is approximated 
by a thermal bremsstrahlung spectrum with a temper- 
ature of ~10-30 keV, depending whether in- or out-of- 
eclipse (Kawai et al. 1989; Cherepashchuk et al. 2005). 
Basically, hydrodynamical radiative plasma flow models 
are found to successfully apply to the data (Brinkmann 
et al. 1991; Kotani et al. 1996). Additional complex fea- 
tures are detected from the XMM spectra, however, which 
could be Compton scattered emission from the jet base 
(Brinkmann et al. 2005) or an iron-K absorption edge due 
to partial covering (Kubota et al. 2007). From the width 
of an eclipse in the 25-50 keV band with INTEGRAL, 
Cherepashchuk et al. (2007) and Krivosheyev et al. (2009) 
propose that a hot extended corona around the accretion 
disk is responsible for the hard X-ray emission via ther- 
mal Comptonization with a temperature of ~20 kcV. High 
quality X-ray spectra covering the broad band are critical 
to establish the interpretation of the high energy spectra 
of SS 433. 

Simultaneous, multi- wavelength observations of this 
source are particularly important to reveal the structure of 
the jets. Previous studies revealed that X-rays are emitted 
from a region close to the base of the jets (~ 10^°~^^ cm; 
Kotani et al. 1996), while the optical emission originates 
far away from the jet base at a distance of ~ lO^"'^^^ cm 
(Fabrika 2004). Radio synchrotron emission is observed 
on a larger scale (> 10^^ cm). Comparison of (quasi) si- 
multaneous data enables us to trace the evolution of the 
same matter traveling along the jets on a spatial scale of 
many orders of magnitudes. 

For this purpose, we organized a multi-wavelength ob- 
servation campaign in the period from 2006 March 29 to 
April 15. The campaign includes two Suzaku observa- 
tions performed in- and out-of eclipse phase and a large 
set of optical spectroscopy, optical/IR photometry and ra- 
dio observations. As for details of this campaign, refer to 
the papers by Kotani et al. (2006, 2008). In this paper, 
we present the results from the Suzaku observations and 
those from optical spectroscopy. Suzaku, the 5th Japanese 
X-ray satellite, covers the wide energy band from 0.2-600 
keV. The large effective area below 10 keV with good en- 



ergy resolution, and high sensitivity above 10 keV, pro- 
vides us the best opportunity to investigate time variabil- 
ity of the wide band spectrum of SS 433. § 2 describes the 
observations and data reduction. The results based on the 
Suzaku data of out-of-eclipse are summarized in § 3, and 
those of in-eclipse are in § 4. We discuss our results in § 5, 
including the comparison of the X-ray and optical data. 

2. Observations and Data Reduction 

2.1. Suzaku 

We observed SS 433 with Suzaku on two occasions, on 
2006 April 4-5 and April 8-9 for a net exposure of Ri40 
ks each. Suzaku (Mitsuda et al. 2007) carries the four 
X-ray Imaging Spectrometers (XIS-0, 1, 2 and 3) coupled 
with the X-ray Telescopes (XRTs) and the Hard X-ray 
Detector (HXD), which covers the 0.2-10 keV and 10- 
600 keV bands, respectively. The XIS-0, 2, and 3 are 
front-side illuminated CCDs (XIS-FI), while the XIS-1 is 
a back-side illuminated one (XIS-BI). The energy reso- 
lution of the XISs is wl30 eV (FWHM) at 6 keV. The 
first observation was performed at the orbital phase of 
= 0.97 (as calculated from the ephemeris by Goranskii 
et al. 1998), during the midst of an eclipse by the compan- 
ion star, while the second was at = 0.26 after the source 
recovered from the eclipse. From our optical observations, 
the precession phase is found to be V' — 0.6 (as calculated 
from the ephemeris by Gies et al. 2002), when we see the 
accretion disk in almost an edge-on view. At this preces- 
sion phase, the directions of the twin jets with respect to 
the line-of-sight are inverted from those observed at most 
of the precession phases, and both show redshifts. In this 
paper, we call the jet with a larger redshift as "red" jet and 
the other as "blue" jet. The XIS was operated in the nor- 
mal clock/editing mode without an window option. The 
target was observed at the XIS nominal position. Table 1 
gives the log of the X-ray observations. 

2.1.1. XIS data 

We analyze the XIS data in a standard manner with the 
HEAsoft^ version 6.6.2 software package. Photon events 
are extracted from a circular region centered at the target 
with a radius of 6 mm (4'. 35). The background is taken 
from the surrounding annulus region. For spectral anal- 
ysis, we make the response matrix file (RMF) and ancil- 
lary response file (ARF) using xisrmfgen and xissimarfgen 
(Ishisaki et al. 2007), respectively. The spectra of all the 
XIS sensors are combined to improve the statistics, to- 
gether with the corresponding responses. We screen out 
events suffering from telemetry saturation, using GTI fil- 
ters provided by the Suzaku team^. 

2.1.2. HXD /PIN data 

The HXD is a non-imaging coUimated X-ray sensor con- 
sisting of two main detection parts; silicon p-intrinsic-n 
diodes (hereafter PIN) and gadolinium silicate crystals 
(Ce-doped Gd2Si05; GSO). They cover the energy range 
of 10 - 70 keV (PIN) and 40 - 600 keV (GSO) with low 

^ http:/ /heasrc. nasa.gov/lheasoft 

^ http: / /www-cr.scphys.kyoto-u.ac.jp/mcnbor/hiroya/xis/gtifilc.html 
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background levels (Takahaslii et al. 2007). In the present 
study, we concentrate on the PIN data, since the signals 
from SS 433 above 40 keV is below the sensitivity of the 
GSO. Since the observations were performed in the nor- 
mal mode of the HXD, we simply extract the spectra from 
the PIN event data, following the standard analysis pro- 
cedure'^. We subtract the non X-ray background (NXB) 
from the data, utilizing the modeled NXB distributed by 
the Suzaku team ^. The version keyword of the NXB 
files is METHOD=LCFITDT (i.e., so-called the "tuned back- 
ground" ) , which gives the best rcproductivity of the PIN 
NXB (0.34% systematic error in the 15-40 keV for a 40 
ksec exposure; Fukazawa et al. 2009). 

In addition to the NXB, two kinds of diffuse/extended 
background should be subtracted before the spectral 
analysis; the cosmic X-ray background (CXB) and the 
Galactic X-ray background (GXB). The GXB is partic- 
ularly important because a target source is close to the 
Galactic plane. Since the surface brightness of the GXB 
varies field to field, we estimate its spectrum using the 
PIN data in a near-by field, {l,b) ~ (28°.5, -0°.2) (OBSID 
500009020), located «11°.5 away from SS 433. This 
region is free from bright X-ray sources and has been 
well studied with Chandra, from which hundreds of faint 
point sources are detected below 10 keV (Ebisawa et al. 
2005; Revnivtsev & Sazonov 2007). They cannot be re- 
solved with the PIN, and hence contribute to the PIN 
signal as a part of the GXB emission. To estimate the 
summed GXB-I-GXB flux, we analyze the PIN spectra in 
this region by subtracting the NXB in the same man- 
ner as above. We find the GXB + CXB spectrum can 
be well fit with a single power-law model in the 12-40 
keV band and obtain the best-fit P = 2.1 and flux of 
2.1 X 10"^^ erg cm~^ s~^. Thus, assuming that the GXB 
(and the CXB) is the same between this region and that 
around SS 433, we subtract the modeled GXB-I-CXB spec- 
trum from the PIN data of SS 433, in addition to the NXB. 
Figure 1 shows the total and source spectra of SS 433 af- 
ter the background subtraction for the eclipse and out-of- 
eclipse phase. The contribution of the NXB, and the sum 
of the GXB and the CXB arc separately plotted. We find 
that the flux level of the GXB-f-CXB is only 0.3 times that 
of the NXB. To evaluate the possible uncertainties in the 
estimated GXB fiux, we refer to the spatial distribution 
of the GXB obtained from the Galactic ridge survey with 
the INTEGRAL/IBIS (Krivonos et al. 2007), which covers 
almost the same energy band as the HXD/PIN. Wc con- 
firm that the systematic errors by using the different field 
separated by 10 degrees are smaller than the statistical 
errors, and hence do not affect our conclusions. Figure 2 
shows the light curves of the HXD/PIN in the 12-40 keV 
from which the NXB is subtracted. 

2.2. Optical Spectroscopy 

We organized quasi-simultaneous optical spectroscopic 
observations of SS 433 using ground-based telescopes 

^ http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
* http:/ /www. astro. isas.ac.jp/suzaku/analysis/hxd/pinnxb/ 
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around the epochs of the Suzaku observations, aiming 
to trace the evolution of the jet structure from the 
X-ray and optical data. From 2006 April 2 to 12, 
four observatories participated in this campaign: the 
Special Astrophysical Observatory (SAO), the Gunma 
Astronomical Observatory (GAO), the Nishi-Harima 
Astronomical Observatory (NHAO), and the Asiago 
Observatory of the University Padova. Table 2 gives the 
logs of the optical observations. Wc observed SS 433 with 
the 6-m Russian telescope^ with the Multi Pupil Fiber 
Spectrograph ^ at the SAO (PI: Fabrika). At the GAO, 
the 150 cm telescope equipped with GLOWS^ (Gunma 
LOW Resolution spectrograph and imager, i?=400-500) 
was used (PI: Kinugasa). At the NHAO, having the 2- 
m Nayuta telescope®, we employed the optical spectro- 
graph (Ozaki & Tokimasa 2005) with a low dispersion 
mode {R ~1000) for our observations (PI: Ozaki). The 
Galileo 122-cm telescope was used at the Padova-Asiago 
Observatory (PI: lijima, i?=:800). Wc also took additional 
spectral data of SS 433 at the GAO about 1-month be- 
fore and 1 and 2-months after the Suzaku observations. 
From these data, we determine the Doppler shift of both 
receding and approaching jets using the Ha emission lines. 
Note that we were not able to detect Ha originating from 
the jets on April 2 (at GAO) and Aprfl 12 (at NHAO) due 
to the limited signal-to-noisc ratio. 

3. Suzaku Results I. (out of Eclipse) 

In this section, wc present the results obtained from 
the Suzaku data on 2006 April 8-9 (out of eclipse phase). 
For spectral fitting, we use the XSPEC package^ version 
11.3.2ag (Arnaud et al. 2006). 

3.1. Light Curves 

Figure 3 shows the summed XIS light curves of the four 
sensors in the three energy bands, 0.5-2 keV, 2-5 keV, and 
5-10 keV for the XIS-FI (0.3-2 keV, 2-5 keV, and 5-8 keV 
for the XIS-BI). As noticed, significant variability is found 
especially in the hard band above the 5 keV. We confirm 
that this tendency also appear in hard X-ray band light 
curve in the 12-40 kcV band obtained with the HXD/PIN 
(see Figure 2). For the subsequent analysis, we divide the 
time region into two, epoch A (the first half), when the 
intensity and spectrum was almost stable, and epoch B, 
when the spectrum exhibited a drastic change. Figure 4 
shows the XIS spectra in epoch A and B. It is seen that 
the average spectrum, in particular the iron-K emission 
line profiles, is very different in the two epochs. 

3.2. Individual Spectral Analysis 

In this subsection, we study the properties of emission 
lines originating from the blue and red jets and their evolu- 
tion based on only the XIS data. To examine the spectral 

^ http:/ /wO.sao.ru/Doc-cn/Telescopcs/bta/dcscrip.htm 

^ http:/ /wO. sao.ru/Doc-cn/TcIcscopcs/bta/instrum/ 

^ http: / /www. astron.pref.gunma.jp/instrumcnts/dcvicc_bcldsi. html 

* http:/ /www. nhao.go.jp/nhao/tclcscopcs/nayuta. html 

^ http:/ /heasarc. nasa.gov/docs/xanadu/xspcc/ 
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variability, we divide the one day data into eight parts 
(data number 1 to 8) with an equal interval, and perform 
spectral fitting to each spectrum. To focus on the iron-K 
band by making the best use of the energy resolution of 
the XIS, we only analyze the 5.0-10 keV band. 

The fitting model consists of a bremsstrahlung contin- 
uum and 11 emission lines, modified by the interstellar 
absorption. It is expressed as a function of energy E: 

^-a{E)Nn ^ [Bremss(fcT) + Fci Ka;j=o 
+(Fexxv Ka + Fexxvi Ka + Nixxvii Ka 
+Fexxv K/3 + Fcxxvi K/S)^^^^^, 
+ (Fe XXV Ka + Fe xxvi Ka + Ni xxvii Ka 
+FexxvK/3 + FexxviK/3)^=^^,„J, 

where cr{E) is the photo-electric cross section and Nn is 
the hydrogen column density toward the source. We fix 
A^H at 1.5 X lO^^cm"^, which is converted with the for- 
mula by Predehl & Schmitt (1995) from Av=8, a typical 
extinction for SS 433 (Fabrika 2004). In this study, for 
simplicity, we adopt a single temperature bremsstrahlung 
model for the continuum, considering the limited energy 
band; although a more complicated model is applied for 
the analysis of the broad band spectra in the next sub- 
section, the results on the emission lines are little af- 
fected. The model includes 11 narrow Gaussian lines, five 
pairs of Doppler-shifted lines from Fexxv Ka, Fexxvi Ka, 
Ni 

XXVII Ka, Fexxv K/3 and Fexxvi K/3 from the two jets, 
and one stationary line of FejKa. The Doppler shifts of 
the red and blue jets are free parameters that are linked 
together among all the lines originating from the same jet. 
The velocity dispersion (Icr width) of the emission lines 
is fixed at 1500 km s~^, which is found by Lopez et al. 
(2006) when the precession phase is the closest to ours 
among available Chandra HETGS observations. 

First, we fit the first four spectra (epoch A) simulta- 
neously, which exhibit little variability. In this fit, the 
intensities of the 6.4 keV FciKa lines are linked one an- 
other among all the spectra. Then, we perform spec- 
tral fitting to the last four spectra (epoch B) separately. 
For these spectra, we fix the intensity of the FejKa line 
at the best-fit value obtained from the epoch A data, 
6.9 X 10~^ photons"^ cm~^, since five emission lines (in- 
cluding four jet lines) are blended around 6.4 keV and it 
is difficult to separate its contribution in epoch B. 

Figure 5 plots the time variability of the center energies 
of the Fexxv Ka lines emitted from the twin jets. As 
noticed, we detect rapid variability of the Doppler shifts 
from data number 6 to 8, with a corresponding rate of 
dz/dt wO. 019/0. 33 day~^. Below we analyze the summed 
spectra only from data number 1-4 (epoch A) , where the 
change of the line energies is smaller than 20 eV. 

3.3. Broad Band Spectra 

We have obtained the simultaneous, broad band X-ray 
spectrum of SS 433 covering up to 40 kcV with the XIS 
and HXD /PIN . This gives us an ideal opportunity to in- 
vestigate the origin of the continuum emission. In this 
subsection, we perform a simultaneous fitting to the XIS 
and PIN spectrum in the 5-10 keV and 12-40 keV band, 



respectively. It is beyond the scope of this paper to ap- 
ply a fully analytic spectral model of the jets that includes 
both continuum and line emission by self-consistent calcu- 
lation. Instead, wc only consider a continuum model, over 
which the 11 emission lines are superposed as independent 
parameters. 

We adopt a multi-temperature bremsstrahlung model 
for the continuum; assuming a power law dependence of 
the temperature of the jet T on the distance from the jet 
base r, we can compute a model spectrum according to the 
expected relation between the differential emission mea- 
sure and temperature. Here wc consider a stationary jet 
model consisting of radiatively cooling, expanding plas- 
mas (e.g., Kotani et al. 1996). The jets are assumed to 
move in a truncated cone with a constant velocity. Hence, 
the density changes as n{r) (xr~^. This model predicts ap- 
proximately T oc r~'^^''~^\ where 7 is the adiabatic index 
(5/3 for non relativistic particles). Since we cannot well 
constrain the power law index from our data, we assume 
T oc r^^^^ according to this model. The free parameters 
are the innermost temperature Tq and its normalization. 

We also consider a reflection component from cold mat- 
ter, which is indicated by the presence of the fluorescence 
6.4 keV iron K emission line, although it has been ne- 
glected in most of previous studies. For this we adopt the 
reflection code by Magdziarz & Zdziarski (1995), avail- 
able as reflect model on XSPEC applicable for any inci- 
dent spectrum. The free parameter is the the reflection 
strength R = n/2Tr, where is the solid angle seen by the 
X-ray emitter (i.e., the jets). We fix solar metal abun- 
dances and an inclination angle of 79°, assuming that the 
reficction mainly takes place in outer parts of the accretion 
disk. 

The fitting results for the epoch A spectra are shown in 
Figures 6 and 7. The best-fit parameters are summarized 
in Table 3. In the analysis, the relative normalization 
between the XIS and HXD/PIN is fixed at 0.923, based 
on the calibration using the Crab Nebula (Ishida et al. 
2007). Wc obtain the innermost jet temperature of Tq = 
27 ± 2 keV. To make the reflection parameter physically 
meaningful, we limit i? < 1 in the fitting and finally obtain 
R= l.Olg 4]^. The equivalent width of the 6.4 keV line with 
respect to the reflection component turns out to be 1.0 
keV at the best-fit value, well consistent with theoretical 
calculations (Matt et al. 1991) within a possible range of 
viewing angle and an iron abundance. 

4. Suzaku Results II. (in Eclipse) 

The same analysis as in the previous section is applied 
for the data of the first Suzaku observation (in-eclipse) . 
Figure 8 shows the X-ray light curves on 2006 April 4- 
5 in the three energy bands. A long term variability is 
clearly seen, as expected from the eclipse of the jets by the 
companion star. Similarly to the out-of-eclipse data, we 
make 8 time-sliced spectra and perform spectral fitting to 
each one, using the same spectral model described in § 3. 
In this period, since the spectra around 6.4 keV arc very 
crowded and hence the intensity of the Fei Ka line can- 
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not be determined, we fix it at 3.7 x lO^^photons^^ cm^^, 
a value that is estimated from the analysis of the broad- 
band spectra with a reflection component, as described 
below. Figure 9 shows the evolution of the line center en- 
ergies of the FexxvKa lines. We do not detect significant 
variability of the Doppler shifts which were observed in 
the out-of-eclipse phase (Figure 5). 

We also make simultaneous fit to the XIS and 
HXD/PIN spectra to constrain the continuum model, us- 
ing the same multi-temperature jet model with a reflection 
component as described in the previous subsection. We 
limit i? < 1 and flx the equivalent width of the fluorescence 
Fei Ka line with respect to the reflection component at 
1.0 keV, based on the result of the out-of-eclipse data, 
since it is difficult to constrain the absolute intensity from 
the data. We obtain Tq = 21 ± 1 keV and R = l.OtEJ ^g. 
Figures 10 and 11 show the data and best- fit model. The 
best fit parameters are summarized in Table 3. 

5. Discussion 

5.1. X-ray and Optical Doppler Shift Curves 

Figure 12 shows the Doppler shifts of the SS 433 jets 
measured with the Suzaku XIS (circles) and optical spec- 
troscopy (crosses). The two dot lines correspond to the 
162.15 days period sinusoidal curves for the blue and 
red jets expected from the precession motion (Gies et al. 
2002), which is shifted by -5.5 days to fit the observational 
points. From this figure, it is confirmed that the Suzaku 
observations were performed at a precession phase when 
the jets axis are almost perpendicular to the line of sight. 

Figure 13 is a blow-up of Figure 12 around the Suzaku 
observation epoch. It is clearly seen that the observed 
Doppler shift curves both from the X-ray and optical data 
show strong deviation from the precession curve with a 
162.15 days period. Except for the data around 2006 April 
8 (MJD 53834) , both amplitude and periodicity of this de- 
viation are consistent with the nodding motion of the jets, 
expressed as a sinusoidal curve with a 6.28 days period and 
a semi-amplitude of Az w 0.012 at our precession phase 
(Fabrika 2004). 

As mentioned above, in the Suzaku XIS data on 2006 
April 8, we detect rapid variability of the Doppler shifts, 
dz/dt « 0.019/0.33 day^^. This is larger than those ex- 
pected from the precession and/or nodding motion of the 
jets; the maximum variability from the nodding motion 
is dz/dt wO.004/0.33 day-^ (Fabrika 2004). Hence, we 
interpret that our Suzaku result corresponds to a jitter 
which was previously reported only from the optical spec- 
troscopy (e.g., lijima 1993), observed for the first time in 
X-rays. The origin of these phenomena will be discussed 
in the next subsection. 

Utilizing our quasi simultaneous X-ray and optical spec- 
troscopic data, we can determine a time lag between the 
epochs when same jet material emits X-rays and opti- 
cal lights, assuming that the Doppler shift remains con- 
stant in traveling along the jet. Based on the optical data 
points on MJD 53834.8 and previous X-ray curves, the 
lag in the red and blue jet is estimated to be 0.64±0.08 



day and 0.38±0.08 day, respectively. Assuming that the 
jets speed is constant at 0.26c, we estimate the distance 
of the optical jets measured from the X-ray jets to be 
/opt = (4.3±0.6) X 10^4 cm (red) and /opt = (2.6±0.6) x lO^^ 
cm (blue). Thus, the Ha emitting region may not be at 
the symmetric location between the twin jets, depending 
on the efficiency of the cooling determined by the plasma 
density. The order of the jet size is well consistent with 
previous estimates (Fabrika 2004). 

The jitter observed in X-rays on 2006 April 8-9 is con- 
firmed in the optical spectra taken after this event. Since 
the length of the optical jets is 1000 times larger than that 
of the X-ray jets, fast time variability observed in X-rays 
is significantly smeared out in the optical band. Thus, 
soon after a rapid change of Doppler shifts is observed in 
X-rays, the optical emission lines become broad (if un- 
resolved) , by including the contribution from jet material 
ejected at different epochs. This is called a "projection ef- 
fect" , previously observed from nodding motions (Borisov 
& Fabrika 1987). Thus, we investigate the change of the 
width of Ha lines originating from the jets before and after 
the second Suzaku observation (out-of-eclipse). The opti- 
cal spectra observed at the GAO are plotted in Figure 14. 
Wc find that while the Ha line width of the blue and red 
jets was 40 A (FWHM) and 58 A on Aprfl 7, the fines 
became significantly broader with a FWHM of 188 A and 
315 A, respectively, on April 9. This change of the Ha line 
width corresponds to Az=0.023 (blue jet) and Az=0.039 
(red jet), which is larger only by A2;=0.005 than the to- 
tal observed change of the Doppler shifts of X-ray lines. 
Thus, the broadening of the Ha lines seen in the April 9 
spectrum is mainly attributable to the rapid change of the 
Doppler shifts occurring close to the jet base. 

5.2. Origin of the X-ray Jitter 

We have detected the fastest variability of Doppler 
shifts of X-ray emission lines ever observed from SS 433, 
which can be interpreted as jitter motions of the jets, not 
the nodding motions where much slower change of red- 
shifts is expected. Assuming that the twin jets are always 
symmetric, we can estimate the orientation and intrinsic 
speed of the jets from the observed Doppler shifts of blue 
and red components. The results are plotted in Figure 15. 
As noticed from the figure, we detect a significant change 
in both angle and speed of the jets during the latter inter- 
val. From the former half toward the end of the observa- 
tion, the jet speed increased by a factor of (8 ±2)%. Recall 
that in the same epoch, the X-ray spectra became softer 
and hence the temperature of the jet decreased. This anti- 
correlation between the jet speed and temperature must 
be explained in theoretical models of jets. 

Considering the smooth evolution of flux and line center 
energy, we infer that the X-ray jitter originates from a con- 
tinuous change of the jet kinetic parameters, not reflect- 
ing the production of discrete jets with largely different 
Doppler shifts. Since the length of the X-ray jets is ^ 10^^ 
cm corresponding to a traveling time of only ^ 100 sec, 
we simply measure the flux-weighted average of Doppler 
shifts of jets emitted during the integration time (~ lO'* 
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sec) of each spectrum. If jets with very different Doppler 
shifts are emitted within each exposure, we would expect 
broadening of the emission lines. We check the line widths 
of the Fc XXV Ka emission lines from the 8 time-sliced 
spectra on 2006 April 8-9. The same fitting model as in 
§ 3 is used, by making the width of Gaussian lines free. 
We find that the width is almost constant over the obser- 
vation at 34 ± 6 eV (the average from the blue jet) and 
36 ±7 eV (red jet), which are comparable to the observed 
line width resolved with Chandra/HETGS (Namiki ct al. 
2003; Lopez et al. 2006). No significant line broadening 
is observed in epoch B. Thus, we favor the continuously 
changing jet model, although we cannot exclude the pos- 
sibility of many discrete jets with slightly different kinetic 
parameters. In the optical jets, by contrast, the evolution 
of emission lines from discrete plasmoids were sometimes 
observed on a time scale of '^0.1-0.3 days (Panferov & 
Fabrika 1993). The relation of these phenomena to the 
X-ray jitter is not clear. 

5.3. Broad Band X-ray Spectra 

We obtain for the first time simultaneous broad band 
spectra of SS 433 covering up to ^40 keV including those 
observed with CCD resolution below 10 keV with Suzaku. 
This gives us a unique opportunity to investigate the ori- 
gin of the spectra not only from the emission lines but also 
from the continuum shape. We find that the spectrum in 
eclipse is consistent with the emission from optically thin 
plasmas in the twin jets as reported in previous works 
(Brinkmann et al. 1991; Kotani et al. 1996), although the 
out-of-eclipse spectrum may need an additional hard com- 
ponent to self-consistently explain both emission lines and 
continuum. 

From the line intensities, we can estimate the innermost 
temperature of the jets based on the radiatively cooling 
jet model, as done in Kotani et al. (1996). The inner- 
most temperature derived from the ratio of the Fexxv 
and Fexxvi Ka emission lines are summarized in Table 3 
as "Jet Temperature" Tq™°. The averaged temperature 
of the blue and red jets is found to be ~13 keV and ~16 
keV, respectively. Due to the confusion with the 6.4 keV 
line under the limited energy resolution, these could be a 
systematic uncertainty in the line intensity of Fexxv- We 
conservatively estimate the possible systematic errors in 
the fine ratio of Fexxv and Fexxvi to be «10%, which 
corresponds to «15% (2-3 keV) in the innermost tem- 
perature. If we take into account this uncertainty, the 
averaged innermost temperature is consistent with each 
other between the out-of-eclipse and in-eclipse phase. In 
our precession phase (edge-on), the region size obscured 
by the outer rim of the accretion disk and that by the 
companion star is similar, and hence there is no large dif- 
ference in Tq'"° between the two Suzaku observations. 

Independently, we are able to constrain the innermost 
temperature Tq- listed in Table 3, directly from the con- 
tinuum shape based on the same jet model, thanks to 
the high quality hard X-ray spectra obtained with the 
HXD/PIN. We find that in the eclipse phase, this Tq value 
(21±1 keV) is close to r^'"'= (16±3 keV, the average of the 



blue and red jets). In the out-of-eclipse phase, however. 
To determined from the continuum (27±2 keV) is signifi- 
cantly higher than that determined from the Fe line ratio 
(13±2 keV). These facts suggest that the out-of-eclipse 
spectrum may be contaminated by an additional hard X- 
ray component that is hidden in eclipse. The spectrum 
must be very hard and/or heavily absorbed since its con- 
tribution should not be large below ~7 keV, where the ra- 
tio of the continuum intensity between the out-of-eclipse 
and in-eclipse phases is similar to that of the iron emis- 
sion lines originating from the jets. The presence of such 
a hard component not directly related to the jets was sug- 
gested by the comparison of in- and out-of-eclipse spectra 
observed with Ginga (Brinkmann ct al. 1991; Yuan et al. 
1995). This may be related to the extended corona around 
the accretion disk discussed in Krivosheyev et al. (2009) or 
a reflection feature of the emission from the inner disk by 
the outer wall of the funnel (Medvedev & Fabrika 2009). 
To constrain its origin, analysis of broad band data uti- 
lizing a fully self-consistent model for the thermal jets 
(including both continuum and emission lines) would be 
required, which we leave for future studies. 
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Fig. 1. The PIN spectra of SS 433 on 2006 April 8-9 (upper, 
out-of-eclipse) and April 4-5 (lower, in-eclipse). From upper 
to lower, the total spectrum, the NXB, the source spectrum 
after background subtraction, and the GXB+CXB. 
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Fig. 2. The PIN light curves of SS 433 in the 12-40 keV 
band on 2006 April 8-9 (upper, out-of-eclipse) and April 4—5 
(lower, in-cclipse). The NXB is subtracted. 
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Fig. 3. The XIS light curves of SS 433 on 2006 April 8-9 
(out of eclipse) in three bands, 0.5-2 keV (upper), 2-5 keV 
(middle), and 5—10 keV (lower). The bin size is 180 sec. The 
arrows define epochs A and B for the spectral analysis. 
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Fig. 4. The 5.0-10 keV spectra of SS 433 in epoch A (upper) 
and B (bottom). The four XISs are summed. See Figure 3 
for the definition of the epochs. 
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Fig. 6. The XIS spectrum in the 5-10 keV band on 2006 
April 8-9 (data number 1-4 only; out of eclipse). The best-fit 
model is over-plotted with separate contribution of the re- 
flection component (lower curve, green) and Gaussians (blue: 
those from the blue jet, red: those from the red jet, green: 
the fluorescence iron-K line at 6.4 keV). The fitting residu- 
als in units of x a-fc plotted in the lower panel. The best fit 
parameters are given in Table 3. 
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Fig. 7. The simultaneous fit to the XIS and HXD/PIN spec- 
tra of 2006 April 8-9 (out-of-eclipse), using a multi-temper- 
ature jet model for the continuum (see text). The best-fit 
models are over-plotted with separate contribution from the 
reflection component and fluorescence iron-K line (green). 
Residuals in units of x '^^'^ shown in the lower panel. See 
Table 3 for the best-fit parameters of the continuum. 



Fig. 5. Time variability of the line center energy of 
FexxvKo determined from the Suzaku spectra on 2008 April 
8 (out-of-eclipse). The horizontal lines show the average val- 
ues determined from the time averaged spectra. 



Suzaku and Optical Spectroscopic Observations of SS 433 



9 



Bin time: 180.0 



> e4 



d o 



p "1 
iri d 
o 



I'/* If* % w M'l 



H')^ #1 «' * 1 



, , , J J ' ' 

I'A i« Vi ,,„s, 



1 2 3 4 5 6 7 8 

— >< — >< — >< — >< — >< — >< >< >■ 

,*J v/r II/. I rf ■••f 



l-.!M <w 



St"* 



2x10* 4x10" 6x10* 8x10* 
Seconds since MJD=53829.5 



10= 



Fig. 8. The XIS light curves of SS 433 on 2006 April 4-5 (in 
eclipse) in three bands, 0.5—2 keV (upper), 2—5 keV (middle), 
and 5-10 keV (lower). The bin size is 180 sec. 
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Fig. 10. The XIS spectrum in the 5-10 keV band on 2006 
April 4-5 (in eclipse). The best-fit model is over- plotted 
with separate contribution of the reflection component (lower 
curve, green) and Gaussians (blue: those from the blue jet, 
red: those from the red jet, green: the fluorescence iron-K line 
at 6.4 keV). The fitting residuals in units of x a-re plotted in 
the lower panel. The best fit parameters are given in Table 3. 
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Fig. 9. Time variability of the line center energy of 
FexxyKct determined from the Suzaku spectra on 2008 April 
4—5 (in eclipse). The horizontal lines show the average values 
determined from the time averaged spectra. 



Fig. 11. The simultaneous fit to the XIS and HXD/PIN 
spectra of 2006 April 4-5 (in eclipse), using a multi-temper- 
ature jet model for the continuum (see text). The best-fit 
models are over-plotted with separate contribution from the 
reflection component and fluorescence iron-K line (green). 
Residuals in units of x ^^e shown in the lower panel. See 
Table 3 for the best-fit parameters of the continuum. 
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Fig. 12. A summary of Dopplcr shifts of the jots observed 
around the observation campaign epoch. The results deter- 
mined from the optical (Ha) and X-ray (FexxvKo) are shown 
by the crosses and circles, respectively. The doted lines cor- 
respond to the 162.15-days sinusoidal curves due to the pre- 
cession motions. 






Fig. 15. The change of the jet parameters observed on 2006 
April 8—9, in terms of the cosine of angle between the 
line-of-sight and jet axis (upper) and intrinsic velocity (lower). 
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Fig. 13. A blow up of Figure 12 around the epoch of the 
Suzaku observations. 
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Fig. 14. The optical spectra of SS 433 observed at the GAO 
on 2006 April 7 (solid line) and April 9 (dashed line). The 
vertical lines denote the line center positions of Ha origi- 
nating from the red and blue jets. There are two Hei lines 
(A6678, A7065) in the figure. 
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Start 


End 




Exposure Remark 


(MJD) 




(MJD) 


(ks) 


Observatory: Suzaku 








2006/04/04 14:40 (53829.6108) 


2006/04/05 12:47 


(53830.5326) 


38.68 Eclipse 


2006/04/08 11:04 (53833.4610) 


2006/04/09 10:59 


(53834.4578) 


40.20 



Table 2. Log of Optical Spectroscopic Observations 



Start 



(MJD) 



Exposure 
(s) 



Remark 



Telescope: BTA 6 m. Observatory: SAO. PI: S. Fabrika. 
2006/04/06 00:45:35 (53831.0317) 590 



Telescope: 150 cm. Observatory: Gunma. PI: K. Kinugasa. 
2006/03/07 19:57:47 (53801.8297) 600(the sum of 120x5) 



2006/04/02 18:47:00 
2006/04/05 19:15:28 
2006/04/06 18:39:53 
2006/04/07 18:40:42 
2006/04/09 18:38:27 
2006/05/03 18:16:07 
2006/05/20 17:35:43 



(53827.7826) 
(53830.8024) 
(53831.7777) 
(53832.7783) 
(53834.7783) 
(53858.7612) 
(53875.7441) 



900(the sum of 180x5) 
900(the sum of 180x5) 
900(the sum of 180x5) 
900(the sum of 180x5) 
1080(the sum of 180x6) 
720 (the sum of 180x4) 
900 (the sum of 180x5) 



Line not detected. 



Telescope: Nayuta 2 m. Observatory: Nishi-Harima. 


PI: S. Ozaki. 


2006/04/03 18:04:52 (53828.7534) 


1800 


2006/04/05 17:55:00 (53830.7465) 


1800 x 2 


2006/04/07 17:48:03 (53832.7417) 


592, 540 


2006/04/12 08:07:16 (53837.3384) 


1800 Line not detected. 


Telescope: 122 cm. Observatory: Padova-Asiago. PI: 


T. lijima. 


2006/04/12 03:12:57 (53837.1340) 


1200 


2006/04/14 02:37:01 (53839.1090) 


1200 



12 



Kubota et al. 



Table 3. The best fit parameters of tfie eontinuum emission 
obtained by a simultaneous fit to tlic XIS+HXD spectra. Tlic 
uneertaintios refer to statistical errors at 90% confidence lim- 
its for a single parameter (but Icr for Jet Temperature Tq'"°). 



Parameter April 8-9 April 4 -5 
Continuum 

To (keV) 27.0tl l Sl.lij i 

R* 1.00_o.4i 1.00_o.i3 

EW** ^ i-OO-o.ll 1.00 (fixed) 

Line Flux: 10~''photons^^cm~'^ 

(rest-frame line energy) Red Blue Red Blue 

FexxvKa (6.698 keV) 40.5^^-^ 45.8;^'^ 17.90l:J]-^^ 19.6lt^-^^ 

FcxxviKa (6.965kcV) 13.9l^-g 12.8^0^5 5.96 ±0J1 S.lllJ^^^ 

Ni XXVII Ka(7.798keV) lO-Sl}'^ U.StH 4.95^°'^^ 5-04l^'^^ 

FcxxvK/3 (7.897keV) 4.7^^-^ 3.8l?-^ 2.8l?f 2.27l°-^^ 

FexxviK/3 (8.210keV) l-7ti-8 3.64 ±0.94 1.90t|^ji^ ^-^^-otl 

FciKa (6.399keV) ' 6.9^2 7 3.7 (fixed) 



Jet Temperature*** T^''^" (keV) 


i4.5±i:^ i2.23±[;:r. 


U.Uto.fs 17.10; 


I-U.93 
-0.87 


Rcdshift 


o.o66i±;J:r^ o.oo64±;i:r^ 


0.0559i;j[;^o^ 0.0318^ 


-0.0U06 
-0.0005 


xVd.o.f 


810/749 


704/698 





*); R = Q/2iT, where Q is the solid angle of the reflector seen from 
the emitter (i.e., jets). 

**): The equivalent width of an Fcj Ko line with respect to the 
reflection component. 

***): Determined from the line intensity ratio between FexxvKa 
and FexxviKa based on a multi-temperature jet model (see text). 



